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Diffuse far-infrared and ultraviolet emission in the 
NGC4435/4438 system: tidal stream or Galactic cirrus? 
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ABSTRACT 

We report the discovery of diffuse far-infrared and far-ultraviolet emission projected 
near the interacting pair NGC4435/4438, in the Virgo cluster. This feature spatially 
coincides with a well known low surface-brightness optical plume, usually interpreted 
as tidal debris. If extragalactic, this stream would represent not only one of the clearest 
examples of intracluster dust, but also a rare case of intracluster molecular hydrogen 
and large-scale intracluster star formation. However, the ultraviolet, far-infrared, Hi 
and CO emission as well as the dynamics of this feature are extremely unusual for tidal 
streams but are typical of Galactic cirrus clouds. In support to the cirrus scenario, we 
show that a strong spatial correlation between far-infrared and far-ultraviolet cirrus 
emission is observed across the center of the Virgo cluster, over a scale of several 
degrees. This study demonstrates how dramatic Galactic cirrus contamination can be, 
even at optical and ultraviolet wavelengths and at high galactic latitudes. If ignored, 
the presence of diffuse light scattered by Galactic dust clouds could significantly bias 
our interpretation of low surface-brightness features and diffuse light observed around 
galaxies and in clusters of galaxies. 
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1 INTRODUCTION 

The study of the diffuse low surface-brightness light around 
galaxies and within clusters represents a unique tool for un- 
derstanding the assembly history of structure in the uni- 
verse. The detection of faint st ellar tails associated wit h 
apparently undisturbed systems ([McConnachie et al.ll2009h , 
as well as the pre sence of diffuse light in galaxy clusters 
(Mih os et al J 12005), are among the clearest pieces of evi- 
dence supporting a hierarchical picture of galaxy formation. 
In the past, the long integration times (and accurate flat- 
fielding) required to detect faint features have hampered 
the search for diffuse extragalactic light. In the next few 
years, with the advent of dedicated deep optical photometric 
surveys such as Pan- Starrs and the Next Generation Virgo 
Cluster Survey, will we eventually obtain an unprecedented 
view of the complex low surface-brightness universe, reach- 
ing very high sensitivities (/i^28 mag arcsec -2 ). The main 
challenge for future investigations will likely be discrimi- 
nating between extragalactic emission and scattered light 
from Galactic cirrus. At low surface-brightness (/is ^27 mag 
arcsec -2 ), scattered light from cirrus dust clouds is clearly 
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detected at optical wavelengths even at high galactic lati- 
tude (e.g. jGuhathakurta &; TvsoiJ ll989). This might be er- 
roneously interpreted as tidal debris or intracluster light. To 
show how difficult it is to discriminate between extragalactic 
and Galactic low surface brightness features, here we present 
the properties of a plume projected near the NGC4435/4438 
pair in the Virgo cluster. This plum e is usually interpreted 
as a tidal stream ([Mihos et al.ll2005h . 

NGC4438 (Arp 120) is the most dramatic example of 
a tidally disturbed galaxy in the Virgo cluster. Dynam- 
ical models suggest that NGC4438 has been perturbed 
by a high-velocity (^800 km s -1 ) interaction with the 
companion SO NGC4435 (e.g., IVollmer et al.1 12003 ). How- 
ever, recent Ha observations favour a different scenario 
in which NG C4438 has collided with the giant elliptical 
galax y M86 (|Kennev et al.l l2008h . Deep optical investiga- 
tions ([Mihos et al. 20051 ) also reveal the presence of a faint 
stellar plume apparently extending from NGC4435. This fea- 
ture has been interpreted in the past as a tidal debris, but 
it is not clear how it would fit into the NGC4438-M86 in- 
teraction scenario. In this Letter, we show that the multi- 
wavelength properties of the optical plume are more similar 
to Galactic cirrus clouds than to tidal debris and that this 



© 0000 RAS 



2 L. Cortese et al. 




Figure 1. From left to right, the Spitzer MIPS 160/zra, 24/xra, V-band (Mihos et al. 2005) and GALEX-FUV images are shown. In the 
right panel, the optical (dashed red) and FUV (solid cyan) contours are superposed to the 160/xm. The contour levels are 25.5, 26.5, 27.5 
mag arcsec -2 and 29.8 AB mag arcsec -2 in V and FUV, respectively. The 24/im and FUV images have been smoothed using a gaussian 
filter with a =3 pixels in order to make the plume more visible. 



feature might not be associated with the NGC4435/4438 
system at all. 



2 THE DATA 

Spitzer observations. We obtained Spitzer MIPS 24, 70 
and 160 fim observations of the NGC4435/4438 system from 
the Spitzer Science Archive. These data are all scan maps 
taken as part of a single proposal (P30945). The image 
frames were created using the MIPS Data Analysis Tools 
(jGordon et all 2005) version 3.10 and using the same tech- 
nique described in lBendo et all <|2009h . The rms in the final 
maps are 0.037, 0.38 and 0.39 MJy sr -1 at 24, 70 and 160 
jum, respectively. The full widths at half maximum (FWHM) 
of the PSFs are 6, 18 and 38 arcsec at 24, 70 and 160 /xm, 
respectively. 

Ultraviolet imaging. The Galaxy Evolution Explorer 
(GALEX) has observed the NGC4435/4438 system in both 
near- (NUV, A = 2316 A, AA = 1069 A) and far- ultraviolet 
(FUV, A = 1539 A, AA = 442 A) bands, as part of the 
Nearby Galaxy Survey. Two different tiles include the two 
galaxies and the stellar tail. We downloaded the final in- 
tensity maps from the MAST archive (GR4/5) and coad- 
ded them, reaching a total exposure time of 7628.35 and 
2967.8 sec in the NUV and FUV bands, respectively. The 
spatial resolution of the final image s is —5". Details abou t 
the GALEX pipeline can be found in lMorrissev et al.1 (2007). 
Single-dish 21 cm Hi line data. Single-dish Hi ob- 
servations of the NGC4435/4438 system were obtained as 
part of the Arecibo Le gacy Fast ALFA (ALFALFA) survey 
([Giovanelli et al.l [2005). The da taset used here comes from 
the first ALFALFA data-release (| Giovanelli et al l2007h , cov- 
ering the center of the Virgo cluster. The total integra- 
tion time is 48 sec beam -1 , providing an rms of —2.2 mJy 
beam -1 at a velocity resolution of — 10 km s -1 . The data are 
gridded into a data cube with V pixels and the spatial reso- 
lution is given by the size of the Arecibo beam, 3.3' x 3.8'. 
Single-dish CO(2-l) observations. To detect CO(2-l) 
associated with the stellar tail, we have carried out a pilot 
observation using the single-pixel A3 receiver (beam-width 
-22") at the James Clerk Maxwell Telescope (JCMTQ). 

1 The James Clerk Maxwell Telescope is operated by The Joint 



We used the full bandwidth of -1.9 GHz (-5OO^U ^19OO 
km s -1 ), thus covering the entire velocity range between 
NGC4435 and NGC4438 (-800 km s -1 ), with a veloc- 
ity resolution of —1.3 km s -1 . A simple ON-OFF observ- 
ing technique was adopted. The ON position roughly cor- 
responds to the peak of 160/xm emission in the stream 
(aj2ooo = 12:27:30.2, £ J2 ooo=+13:12:29), while the OFF po- 
sition was accurately chosen in order to avoid any possi- 
ble contamination from Galactic emission. Data were re- 
duced using the Starlink data reduction packages (KAPPA 
and SMURF). Individual spectra were inspected and then 
coadded to produce a final spectrum with a total on-source 
integration time of 1347s. The rms on the final baseline- 
subtracted spectrum is —0.018° K (expressed in antenna 
temperature). 



3 RESULTS 

In Fig. [1] we compare the Spitzer 24 and 160 /mi images, 
the GAL EX FUV image, an d the deep V-band image ob- 
tained bv lMihos et~aD (|2005h for the NGC 4435/4438 sys- 
tem. They clearly reveal the presence of far- infrared and far- 
ultraviolet diffuse emission associated with the low surface- 
brightness (fjLv —26.5-28 mag arcsec -2 ) optical plume to 
the north /n orth-west of NGC4435 originally discovered by 
iMalinl (119941 ). While the MIPS field of view only includes the 
brightest part of this feature (close to NGC4435), the larger 
area covered by GALEX allows us to trace the whole extent 
of the stream out to a distance of — 13 / from NGC4435. The 
plume has an unusual L-shaped (or triangular) morphology, 
with abrupt and well defined edges to the east and north. 
No resolved sources (e.g., star- forming knots) are visible in 
the plume. 

Interestingly, we first noticed the presence of diffuse 
emission in the Spitzer images, where the intensity of this 
feature is particularly high: —1.5-5 MJy sr -1 and —0.05-0.08 
MJy sr -1 at 160 and 24 jLtm, respectively. A detailed inspec- 
tion of IRAS images reveals that the plume is marginally 

Astronomy Centre on behalf of the Science and Technology Fa- 
cilities Council of the United Kingdom, the Netherlands Organi- 
zation for Scientific Research, and the National Research Council 
of Canada. 
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Figure 2. The Hi contours of the plume superposed on the FUV 
image. The blue contours show the Hi emission detected by AL- 
FALFA (including diffuse Galactic Hi ) across the heliocentric 
velocity range — 9.2< Vq < —19.5 km s _1 . Contours levels are 
7.2 to 8.2 Jy km s — 1 beam — 1 , in steps of 0.1 Jy km s — 1 beam — 1 . 
As discussed in §3, the net Hi intensity of the plume has been ob- 
tained by performing a median subtraction of the Galactic emis- 
sion around the feature. For comparison, the VLA contours of the 
tail published by Hota et al. (2007) are shown in red. We note 
that part of the Hi associated with NGC4438 and IC3355 is not 
visible given the narrow velocity range here investigated. 



visible in the 100 /im image, although the lower spatial 
resolution of IRAS made it impossible to clearly associate 
the plume with the NGC4435/4438 system. No emission 
is detected at 70 /im, implying a 3a upper limit of ~1.14 
MJy sr _1 . The GALEX data reveal that the plume is really 
evident in FUV (hfuv ^29 AB mag arcsec -2 ), whereas 
in NUV it is only very marginally detected at a surface- 
brightness level of unuv ^29.2 AB mag arcsec -2 . 

In addition to dust and stars, this feature appears also 
to contain a significant amount of ato mic hydrogen and CO. 
Using the Very Large Array (VLA) . [Hota et all (|2007l l de- 
tected Hi emission associated with the brightest part of the 
plume. Intriguingly, the Hi has a low recessional velocity 
(rsj —9.5 km s -1 , for a velocity resolution of 20.7 km s -1 ) and 
a narrow velocity width (<40 km s -1 ). Unfortunately, the 
northern and eastern edges of the tail were too close to the 
half-power of the VLA primary beam to be detected. Thanks 
to the higher sensitivity, velocity resolution and larger field 
of view of ALFALFA observations, we are now able to fill this 
gap and to characterize in more detail the Hi properties of 
this plume. By inspecting the cube through a velocity range 
-200< V <1000 km s _1 (i.e., the typical interval expected 
if the plume is related to the NGC4435/4438 system), we 
confirm that the only Hi emission clearly associated with 
the optical plume has a heliocentric velocity V© = — 15±4 
km s _1 and a width =5± 7 km s -1 . In Fig. [2] the AL- 
FALFA and I Hot a et al. (|2007l ) Hi contours are superposed 
to the FUV image of the NGC4435/4438 system. While the 
two datasets fairly agree in the overlapping region, the AL- 
FALFA observations are able to trace the Hi along the whole 
extent of the plume, showing a good match between the Hi 
and FUV morphology of this feature. To estimate its net Hi 
intensity we performed a median subtraction of the Galactic 
emission around the plume. The net flux density varies in the 
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Figure 3. The single-pointing CO(2-l) (solid line) and integrated 
Hi (dotted line) spectrum of the plume. 

range 0.6-1.2 Jy km s _1 beam _1 . Assuming that the Hi fills 
the Arecibo beam, this can be transl ated into a colum n den- 
sity of -1.5-3xl0 19 cm" 2 (see also iHota et~aH l2Q07h . The 
integrated Hi spectrum of the plumjfl is shown in Fig. [3] 
The peak and absorption features on either side of the de- 
tection are residual artifacts of the sky subtraction based on 
the median values in the selected channels. The recessional 
and velocity widths of the plume are confirmed, and even 
better constrained, by our CO(2-l) observations, as shown 
in Fig. [3] CO(2-l) emission is clearly detected (peak signal- 
to-noise ^20) at a heliocentric velocity V© = — 13.9±0.1 
km s _1 and with a FWHM = 1.5±0.2 km s" 1 . 

In summary, the plume shows emission coming not only 
from stars and atomic hydrogen but also dust and CO, 
emerging as an almost unique feature in the panorama of 
low surface-brightness streams and tails discovered so far in 
the Virgo cluster. 



4 DISCUSSION 

What is the nature of this plume? Is it really a tidal stream 
associated with the NGC4435/4438 system? If extragalactic, 
it would represent an extremely rare example of intracluster 
dust and molecular hydrogen stripped by gravitational in- 
teractions. However, the mult i wavelength properties of this 
feature demand a little bit of attention before one can inter- 
pret it as a tidal feature. This is particularly true since the 
plume is detected in far-infrared, where the contamination 
from foreground Galactic cirrus emission might be signifi- 
cant, if not dominant. In the following, we investigate the 
two possible scenarios (tidal stream and Galactic cirrus) in 
order to unveil the most likely origin of this feature. 

4.1 The tidal scenario 



If at the distance of the Virgo cluster (16.5 Mpc, iMei et al.l 
I2007T ). the plume would have a physical size of ^33 kpc (in 
both right ascension and declination), extending up to ^70 
kpc from NGC4438. Its Hi and CO recessional velocities 
(rsj —14 km s _1 ) are roughly consistent with the Hi veloc- 
ity of NGC4438 (V& =104±2 and W50 =244±5 km s~\ 
iGiovanelli et al.l I2007T ) but significantly different from the 

2 We only considered the main Hi feature visible Fig(2] i.e., 
the one clearly associated with the optical, UV and far-infrared 
plume. 
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Figure 4. The Virgo cluster region as seen in FUV (left) and 100 (right). The region shown in Fig. [T] is indicated by the white 
rectangle. The FUV mosaic have been smoothed using a gaussian filter with a =9 pixels in order to highlight the diffuse cirrus emission. 



redshift of NGC4435 (V -801 km s" 1 ). This apparently 
rules out a gravitational interaction between NGC4438 and 
NGC4435: given their large velocity difference, we should ex- 
pect the stripping material to lie at a velocity intermediate 
between the two systems. 

The plume apparently connects NGC443 8 to the dwarf 
irregu lar galaxy IC3355 (V — —17 km s~\ iGavazzi et ah! 
2003), perhaps suggesting a recent interaction between the 
two systems. This seems sup ported by the prese nce of a Hi 
tail associated with IC3355 ([Chung et al.ll2009h . However, 
contrary to what expected in case of a fly-by encounter with 
NGC4438, IC3355's tail points to the west, thus in the op- 
posite direction with respect to NGC4438. Moreover, it is 
still unclear w hether or not IC3 355 is really a member of 
the M86 cloud (jChung et alJl2009h . Thus, although we can- 
not exclude that IC3355 is playing a role, the current data 
do not allow us to further investigate this hypothesis. 

A more plausible scenario is that the plume has been 
created during the collision between NGC4 438 and M86 
V — —244 km s _1 ) as recently unveiled by iKennev et al.l 



(2008). A complex structure of Ha filaments connects the 
two galaxies and spans a velocity range (— 240 < V <+70 
km s _1 ) consistent with the value observed in the plume. 
Moreover, the morphology of our feature (i.e., well defined 
edges on the east with less clear boundaries on the west 
side) would be consistent with the direction of motion of 
NGC4438, from w est to east. The tim e-scale of the interac- 
tion is —100 Myr (|Kennev et al.ll2Q08h . significantly shorter 
than the dust sputterin g time-scale expecte d at the position 
of NGC4438 (-10 8 vn lPopescu et alJfeOOd ). Thus, although 
extremely unusual, the presence of dust in the tail does not 
rule out the tidal scenario. 

However, not all the properties of the plume are entirely 
consistent with a gravitation interaction. Above all, its ve- 
locity width (—2-5 km s _1 ) is significantly smaller than the 
typical values (>20 km s _1 ) observed in tidal tails, candi- 
date tidal-dwarf galaxies, isolated Hi clouds and ram pres- 



sure strijDj^edj^ IKent et aD 

120071 : lOosterloo & van Gorkoml 12005). Although we cannot 
exclude that the cloud is seen 'face-on', we consider this to be 
unlikely. The presence of a clear velocit y gradient in the Ha 
filaments between M86 and NGC4438 (jKennev et alJl2008h 
suggests in fact that the tidal acceleration is not parallel to 
the plane of sky. Moreover, it is difficult to imagine how a 
violent process such as a tidal interaction would be able to 
strip material along a well defined direction and to keep the 
velocity dispersion of the Hi cloud so small. It is therefore 
unclear how tidal debris can maintain such small velocity 
spread on a scale of a few tens of kpc. 

Equally unusual are the UV properties of the tail. The 
FUV emission is in fact diffuse and not patchy with star- 
forming knots as observed in s tar-forming disks, t idal tails 
and extended UV disks (e.g., iThilker et aD l2007h . Extra- 
planar diffuse UV emission is often interpreted as being 
due to light scattered by diffuse dust, in particular when 
associated with th e halos of starburst galaxies (like M82; 
iHoopes et al.ll2005h . Although the Spitzer images reveal that 
dust is surely present in the plume, the lack of a power- 
ful central starburst and the large distance of the tail from 
any strong radiation source apparently exclude the scatter- 
ing scenario. Thus, if extragalactic, the FUV emission must 
come from young stars in the plume. The absence of star- 
forming knots and Ha emission (jKennev et al.ll2008l ) implies 
an age between —10 Myr (i.e., the lifetime of OB associa- 
tions emitting in Ha) and —25 Myr (i.e., the lifetime of B 
stars, supposed t o be responsible for the diffuse UV light 
in galactic disks; IPellerin et al.1 [20071 ) , significantly shorter 
than the time-scale of the interaction but consistent with the 
ag e of the UV tidal ta il to the west of NGC4438 discovered 
bv lBoselli et al.l (|2005h . This might suggest that the plume 
has just entered a post-starburst phase, consistent with the 
fact that its colour is significantly bluer t han the observed 
colou r of NGC4438 (FUV - V -5.2 mag; iGil de Paz et~ai1 
l2007h . However, if this is really the case, it remains unclear 
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to us why the plume stopped forming stars. The very low 
velocity dispersion, the large amount of Hi and the presence 
of molecular hydrogen should allow the cloud to continue 
forming stars. Moreover, the fact that CO is usually asso- 
ciated with dense st ar-forming regions, e ven in tidal tails 
and stripped clouds (Lisen feld et alj[2004 h makes the post- 
starburst scenario unlikely. 

4.2 The Galactic cirrus hypothesis 

A completely different possibility is that the plume is in re- 
ality a Galactic cirrus cloud of gas and dust just superposed 
by chance near NGC4435. In this case, the diffuse FUV and 
optical emission would come from light scattered by the dust 
grains in the cloud. 

Before testing the cirrus hypothesis in detail, it is worth- 
while to examine the probability of a cirrus cloud appear- 
ing superimposed in front of the Virgo cluster. Despite 
its high Galactic latitude (b ^75°), the center of Virgo 
lies in a region that is significantly contaminated by cir- 
rus emission. A ring-like cirrus struct ure centered near M87 
and extending over several degrees (|Brosch et al.l 1 19991 ) is 
clearly visible in the IRAS 100/im images (see Fig.[4l right). 
Since scattered light is often assoc iated with cirrus clouds 
(e.g. JGuhathakurta fe Tvsonlfl989h . we should expect to see 
the same ring at both ultraviolet and optical wavelengths, 
though only at very low surface-brightness (/is ^27-28 mag 
arcsec -2 ). While deep optical data of the whole Virgo re- 
gion are not available, the deep and wide GALEX coverage 
of the Virgo cluster makes it possible to test this hypothesis. 
We thus combined all the available GALEX pointings near 
the center of Virgo having FUV exposure time longer than 
800 sec. The mosaic is shown in the left panel of Fig. [H and 
it likely represents the highest resolution large-scale map of 
FUV cirrus emission made to date. It clearly illustrates the 
presence of diffuse FUV emission tracing the distribution of 
far-infrared Galactic cirrus emission. In addition the high 
spatial resolution of GALEX allow us to see for the first 
time the detailed structure of the cirrus clouds. In order 
to quantify the correlation between FUV and FIR diffuse 
emis sion, we first removed all sources detected by Sextrac- 
tor ([Bertin &; Arnoutsf l996) and then smoot hed the UV im- 
ages t o the resolution of the 'point-source free' Schl egel et al.l 
(1998) IRAS map (~6 arcmin) by applying a Gaussian filter. 
The Pearson correlation coefficient between the FIR and UV 
intensity in each pixel (142x142 arcsec 2 ) is ~0.82, consistent 
with the range of values obtained by previous analysis of cir- 
rus c louds (e.g., Haik ala et al.l [l995l ; ISasseen &; Deharvend 
1996). Thus, there is a reasonably high probability of detect- 
ing diffuse scattered far-ultraviolet (and perhaps also opti- 
cal) light associated with cirrus emission in the direction of 
the Virgo cluster. 

More direct support to the cirrus hypothesis is provided 
by the properties of the plume. Its negative recessional veloc- 
ity, narrow velocity width and Hi column density are com- 
pletely consistent wi th what is observed in Galactic low- 
velocity clouds (e.g., [Stanimirovic et al.| [20Q6^) . In addition, 
the presence of CO emission (and thus of mole cular hydro- 
gen) is commonly observed in cirrus clouds ([Reach et al.l 
1 19981 ). 

Several studies have shown that there is a good correla- 
tion between the 100/xm far-infrared intensity and the Hi col- 



umn density of cirrus cl ouds: e.g., IjlOO/im) /N (HI) ^0.5-3 
MJy sr -1 / 10 20 cm 2 (|Boulanger et alJll996h . At 160/xm, 
the far- infrared to Hi ratio of the plume is ~10 MJy sr -1 
/ 10 20 cm 2 . Assuming a typical 160 to 100 /xm flux ratio ~ 
1.5-3 for cirrus clouds (|Bot et alJ [2009), our feature appears 
to have a significant far- infrared excess. This is however not 
inconsistent with the cirrus hypothesis since Galactic 'far- 
infrared excess clouds' are not rare at high Galactic latitudes 
and, as observed in this case , they are usually associated 
with significant CO emission ([Reach et al.lll998r ). Thus, the 
far-infrared-to-Hl ratio does not allow us to conclusively dis- 
criminate between the two proposed scenarios. This is also 
because our knowledge of the gas and dust content of cirrus 
is mainly based on studies of physical scales of tens of ar- 
cminutes, and it is not yet clear whether we can extrapolate 
them to the physical scale of the plume. 



5 CONCLUSION 

In this Letter, we have investigated the multiwave- 
length properties of the optical plume projected near the 
NGC4435/4438 system. It is very tempting to interpret this 
feature as tidal debris in the Virgo cluster, as this would rep- 
resent an extraordinary case of intracluster dust and molec- 
ular hydrogen as well as large-scale intracluster star forma- 
tion. However, our analysis reveals that the dynamics, mor- 
phology, gas and dust content of this stream are more con- 
sistent with what observed in Galactic cirrus clouds than in 
tidal features associated with cluster galaxies. Whatever its 
real origin, this stream is a fascinating object and it clearly 
highlights how difficult is to discriminate between Galactic 
clouds and extragalactic diffuse light. This will be a great 
challenge for both far-infrared Herschel surveys and deep 
optical investigations of the extragalactic sky. Only a de- 
tailed characterization of the FUV-to-far-infrared properties 
of cirrus at all physical scales might eventually allow us to 
disentangle between the two different scenarios. 
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